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Schwinger-Keldysh canonical formalism for electronic Raman scattering 
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Inelastic low-energy Raman and high-energy X-ray scatterings have made great progress in instru¬ 
mentation to investigate the strong electronic correlations in matter. However, theoretical study of 
the relevant scattering spectrum is still a challenge. In this article, we present a Schwinger-Keldysh 
canonical perturbation formalism for the electronic Raman scattering, where all the resonant, non¬ 
resonant and mixed responses are considered uniformly. We show how to use this formalism to 
evaluate the cross section of the electronic Raman scattering off an one-band superconductor. All 
the two-photon scattering processes from electrons, the non-resonant charge density response, the 
elastic Rayleigh scattering, the fluorescence, the intrinsic energy-shift Raman scattering and the 
mixed response, are included. In the mean-field superconducting state. Cooper pairs contribute 
only to the non-resonant response. All the other responses are dominated by the single-particle 
excitations and are strongly suppressed due to the opening of the superconducting gap. Our formal¬ 
ism for the electronic Raman scattering can be easily extended to study the high-energy resonant 
inelastic X-ray scattering. 

PACS numbers: 78.20.Bh, 78.30.-j, 74.25.nd, 61.05.cc 


I. INTRODUCTION 


Inelastic low-energy Raman and high-energy X-ray scatterings have become powerful tools to study the strong 
electronic correlations in matteriir—. While the instrumental technique of the light scattering is in rapid progress, the 
theoretical study of the scattering spectrum is in less development. Two main difficulties suppress the theoretical 
study of the light scattering off the strongly correlated electrons. The first difficulty stems from the complexity 
of the strongly correlated electrons themselves. We now have no well-defined theoretical formalism for the various 
electronic correlations in such as the high-Tc cuprates, iron-based superconductors and heavy fermions etc., where the 
multiple comparable energy scales and different degrees of freedom are strongly correlated. The other difficulty lies in 
the description of the inelastic light scattering processes. Unlike the single-particle scattering technique such as the 
angle-resolved photoemission spectroscopy (ARPES), neutron scattering etc., Raman and X-ray scatterings involve 
two-step photon-in photon-out processes. The cross section of the ARPES or neutron scattering is determined by the 
scattering correlation function which can be studied in perturbation formalism by the fluctuation-dissipation theorem. 
However, a simple extension of this formalism for the two-photon scattering is fail because the fluctuation-dissipation 
theorem is now invalid. We thus have no reliable perturbation formalism to study the scattering correlation function 
in Raman and X-ray scatterings. 

In this article, we focus our study on the second difficulty. We show that it can be overcome by introducing the 
Schwinger-Keldysh contour time formalism, which has been well established for non-equilibrium physics^^— . In this 
article, we present a Schwinger-Keldysh perturbation formalism to evaluate the cross section of the electronic Raman 
scattering. The formalism for the high-energy resonant inelastic X-ray scattering (RIXS) can be established in a 
similar procedure. 

Our starting point is the differential cross section of the inelastic light scattering. Consider a two-step photon-in 
photon-out scattering as shown schematically in Fig. [1] The incident photon with momentum pi and polarization 
is absorbed by the electrons of the target matter which then emits photon with momentum p f and polarization e f . 
Suppose the initial state of the electrons is |^i) at time ti and the final state after the scattering is |^/) at time tf. 
The scattering probability of this two-photon process is described by 


r(p/e/;piei) = 

4>i4>f 




2 


( 1 ) 


where S{tf,ti) is the time evolution matrix from an initial state jT/) = Ipie^^i) into a final state I'I'f) = Ip/g/^/), 
and Ei is the energy of the electrons in the initial state. Suppose there are N photons in the initial state [piej. 
Among the N photons there are photons scattered. The conservation of the photons in the 

scattering process shows that 


$i(pj,ei)crAt = N r(p/e/;pjej), 

P/e/ 


( 2 ) 
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where a is the effective scattering cross section, = nc = ^ is the current density (or flux) of the incident 

photons (y is volume of the photon field and c is the light velocity), and At = tf — . Since uif = pfC, we have 

Spy = ( 2 ^ 0 )'-^ I iJj'fdujfdfl where dfl the differential solid angle. The double differential cross section with the initial 
and final photon states Ip^eJ and |p/e/) is given by 




dTtdjjjf 




p5^pjj^r(p,eyp.e.), 


( 3 ) 


where q and v are the transferred momentum and energy frequency, respectively, and are defined by 


q = Pi-pf,i2 = u}i-u}f. (4) 

Formula ([31) shows that the differential cross section is proportional to the scattering probability F. The time difference 
At can be canceled by an additional factor At in F which comes from the energy conservation law. Therefore T/At 
can be taken as a scattering rate. 



FIG. 1 : Schematic light scattering from a target matter. |piei) and |p/e/) denote the incoming photon state with momentum 
Pi and polarization ei and the scattered photon state with momentum p/ and polarization er, respectively. 

Suppose the coupling of the electron and the photon field is Vj. Define the total Hamiltonian of the combined 
system as Jif = H + Hp + Vj with H and Hp the Hamiltonian of the electron and the photon system respectively, S 
matrix is given by 


S{tf,ti) = Tte 


i dtviit) 
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( 5 ) 


where Vi{t) = jg time ordering operator. Separate the interaction Vj into 

Vi of linear to A and V 2 of quadratic to A , where A is the photon vector potential. To lowest-order perturbations, 
only the following two expansions of the S matrix contribute to the scattering probability F, 


r=E 

4>i4'f 


,-PEi 

Z 


+S'2|4'/) 


( 6 ) 


where S '!.2 are defined by 

Si = dtVy), (7) 

^2 = dtidt2Tt[Vi{ti)Vi{t2)]. 


Thus the scattering probability F involves three contributions. 


r = ri + r2 + ri2. 


( 8 ) 
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with 


hi 


r 2 


ri 2 


E 

<t>i4>s 

E 

4'i4>f 

E 


.-PEi 

Z 

.-PEi 

z 

,-PEi 

~z 




-2Re 


(9) 


Fi, r2 and ri2 are the so-called non-resonant, resonant and mixed parts of the scattering probability, respectively. 
ri2 describes quantum interference of the resonant and non-resonant scattering processes. The positive or negative 
ri2 comes from the corresponding constructive or destructive quantum interference. 

Since the states of the incident and the scattered photons are defined definitely, the scattering probability F can 
be reduced into a representation of the pure electron system. Now 51^2 matrices can be re-expressed in similar forms 
to Eq. © where the interactions hi,2 are substituted by the reduced ones Vi,2 without photon field involved (details 
and derivation will be shown in the following section). 

The non-resonant scattering probability Fi is determined by the correlation function as 

Fi = ^ ' dtidt2{Vl{h)V2{t2)), (10) 


where (A) = ^^-Tr 


,-PHa 


and V{t) = *i)Ve U)_ the fluctuation-dissipation theorem, Fi can be 


re-expressed into the standard form: 


2 At 

where x(f) is the frequency Fourier transformation of the time-ordered correlation function x(ti,t2) = i 0 {ti — 
l2)([V|(ti), V2(t2)])- Perturbation theory can then be easily introduced to evaluate Fi. This is a standard formalism 
to study the scattering probability in the single-particle scattering technique such as ARPES, neutron scattering etc. 

Because of the time ordering operator It[Vi(fi)Vi(t2)] in S2 matrix, the fluctuation-dissipation theorem is invalid 
to evaluate the resonant F2 and the mixed F12. In most studies of the Raman or the X-ray scattering spectrum, 
F2 and Fi 2 are evaluated from the Kramers-Heisenberg formulaii^, where the perturbation is badly controlled and 
numerical methods are applied. No reliable perturbation formalism is established even for the weakly interacting 
electron system. The time ordering in S2 matrix is the difficulty we should overcome to establish a perturbation 
formalism to evaluate F2 and F12. 


c. 

< ti 

c 

FIG. 2 : Two-branch contour C for time ordering operator Tc. ti and tf are the initial and final times respectively. C = C+UC- 
with an upper time branch C+ '.ti^tf and a lower time branch C_ : tf —>■ ti. If tt —>■ —oo,tf —>■ -too, the contour C is the 
so-called Schwinger-Keldysh contoup^. 

From the picture of a time evolution, the scattering probability F involves two time evolution processes, forward 
time ordering from the initial state Idli) at ti to the final state IdlF) at tf , and backward anti-time ordering from the 
final state |^'f) back to the initial state Introducing an anti-time ordering evolution matrix S, 

( 12 ) 

we can then define the anti-time ordering Si^2 matrices analog to ^1,2 in Eq. ©• Following the time-and-anti-time 
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evolution picture, r2 can be expressed as 

r2 = ^ — — —('I'/|5'2|4'f)(^f|*52|'I'/) 

= ^ [ m[v/(i2)vl(t'i)]rt[Vi(ti)Vi(t2)]), 

where /[(.(/] = dt'idt2 dtidt2 (this abbreviation will be used in the whole article where the time variable 
without or with prime will follow time evolution or anti-time evolution, respectively). Introduce a time contour C 
which describes both the time and the anti-time evolution process, C = C+ U C_, where t £ C+ evolve as ti —>■ t/ 
and t' £ C- evolve as t/ —>■ as shown schematically in Fig. [21 Then r2 can be re-expressed as 


r2 = i / (T,[vJ(t'2)vJ(<;)Vi(ti)Vi(t2)]). 

Here Tc is the contour time ordering operator defined by 


T,[Aiti)B{t2)] 


A{ti)B{t 2 ), if ti >c < 2 , 
±B{t 2 )A{ti), if ti <c t 2 , 


( 13 ) 

( 14 ) 


where >c and <c are defined according to the position of the contour time arguments, latter or earlier in the time 
contour C, and ± are defined for the bosonic or fermionic operator, respectively. 

From a similar derivation, the non-resonant Fi and the mixed F12 can be re-expressed within the contour-time 
formalism as 


ri = - / (rc[v|(<')V 2 (t)]), 


( 15 ) 


ri 2 


Re 


* / (T,[v|(t')Vi(ti)Vi(t2)]) 


( 16 ) 


Introducing the contour time formalism to evaluate the scattering probability F is our principle to overcome the 
difficulty in theoretical study of the two-photon inelastic light scattering. Formulae ©, (USD, (USD and (USD constitute 
the contour time formalism. When the times are set as ti —>■ — 00 and t/ —t -l-oo, this formalism becomes a Schwinger- 
Keldysh formalism which we will use in realistic calculation. 

Introduction of the Schwinger-Keldysh contour time formalism into the study of the resonant inelastic X-ray scat¬ 
tering can retrospect to 1974 by Nozieres and Abrahams^, which was then followed by Igarashi et al in 2006 ^. In their 
formalism the authors focus on the scattering rate, a time derivative of F with respective to tf. This scattering rate as 
a time derivative leads to the broken equivalence of the time variables and thus the formalism obtained is limited in 
study for its complex. In our formalism, the double differential scattering cross section is related directly to the scat¬ 
tering probability F and thus all the time variables are in equivalent symmetry. A path integral functional formalism 
is provided recently by H. C. Lee^id^. While only the resonant scattering process is studied for the resonant inelas¬ 
tic light scattering measurement^, all the two-photon scattering processes including the resonant, non-resonant and 
mixed responses from magnons in antiferromagnetic insulators are considered uniformlyiS. In the article, we present 
an equivalent canonical formalism to study the electronic Raman scattering, where all of the electronic responses, 
the non-resonant charge density response, the elastic Rayleigh scattering, the fluorescence, the intrinsic energy-shift 
Raman scattering and the mixed response, are included in our theory. 

The article is arranged as following. In Section [H we present the principle to establish a Schwinger-Keldysh contour 
time formalism for the two-photon inelastic light scattering. In Section im we show an example how to use this 
formalism to study the electronic Raman scattering off an one-band superconductor. The scattering cross section 
with contribution from all the two-photon processes is studied in details in mean-field approximation. Summary is 
present in Section m In 13 we provide a preliminary introduction to the non-equilibrium contour time formalism^ 
for those who are not familiar with it. 


II. ELECTRONIC RAMAN SCATTERING OFF AN ONE-BAND SUPERCONDUCTOR 

In this section, we show how to study the electronic Raman scattering off an one-band superconductor with the 
Schwinger-Keldysh formalism we have established in SectionlH The contribution from the resonant, non-resonant and 
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mixed responses to the scattering cross section is evaluated in details. All the responses are considered in mean-field 
approximation in superconducting state as an example. 

It should be noted that in superconducting state some low-energy excitations may play roles in the electronic Raman 
responses, which include the Bardasis-Schrieffer bound statesii^— , the longitudinal and transverse phase modesi^i^— , 
the amplitude Higgs modei^ as well as the orbital excitations in multi-orbital superconductorsi^. The Bardasis- 
Schrieffer bound states are difficult to be resolved in experiments due to the small binding energy, the finite life 
time and/or the weak spectrum weight. The longitudinal phase mode which is important for the gauge invariance is 
modified into high-energy plasma by the Coulomb interaction and becomes low-energy irrelevant in Raman responses. 
Since the main features of the Raman responses are dominated by the gapped Cooper pairs in superconducting 
state with a transverse renormalizationi^r— , in this article we make a mean-field approximation with only pairing 
interaction involved. This mean-field approximation is also suitable for simplicity to show how to calculate the Raman 
responses with the Schwinger-Keldysh formalism we have established. The roles of the relevant low-energy excitations 
in the Raman responses are other important issues to be studied in future. 


A. Scattering probability in contour time formalism 

Consider an one-band electron system with Hamiltonian H = Ht + i?/, where Ht and Hi are the free kinetic and 
the interaction part of the Hamiltonian respectively. Ht is given by 

Ht = ~^^Ujdl^dja, (17) 

ijn 


where dimd\^ are the annihilation and creation operators, respectively, of the electron at site i with spin tr. The 
electron-photon coupling can be obtained by considering the gauge invariance of Ht, which leads to an additional 
phase factor for ty and thus 


Ht{A) 


ija 




(18) 


( R -f-R ■ \ 

—- 1 is defined on bond and the charge of electron is —e. Extend the phase factor into second-order 
of A, the electron-photon coupling Vj = Vi + V 2 can be obtained as 


( 19 ) 

qa 

R 2 = ^n“^(-qi-q2)A“(qi)A/'(q2), 

qiCt/S 

where a,f} = x,y,z. /“(—q) and n“^(—q) are the current and the stress tensor operator which couple linearly or 
quadratically to A respectively, and are defined by 


^ ktr 

(20) 

n“^(-q) = ^X]T“^(k,q)4+q^rfka, 

ktr 


u“(k,q) = 

<5 

T“4k,q) = 

(21) 


Note that the electron-photon coupling in (1191) only involves the electron charge degree of freedom. In a more general 
case, the electron-photon coupling should also involve the electron magnetic orbital and spin degrees of freedom^. 
A similar derivation can be done for the multi-orbital electron system such as Fe-based superconductors, where the 
orbital fluctuations may have unusual contribution to the Raman scattering. 
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Let us now derive the scattering probability. Introduce the second quantization of the vector potential A as^i 

' h 


A(q,t) = 


2eQUJc^xV 


eA(q) (aqA(i) + aqAW) > 


where ac^\{t) = aqAe“®‘^‘i^‘ and eA(q) is the polarization vector with A = 1,2 (the polarization can be linear or 
circular). In the below we will assume that the photon energy is polarization independent, i.e., WqA = Wq. When the 
photon states in Eq. dZl)-® are traced out, the non-resonant scattering probability Ei follows as 


Li 



4>i4>i 





2 


where the operator V 2 {t) is defined by 

kcr 


( 22 ) 


(23) 


Here ro(t) and A(k, q) are given by 

ToW . e-“,A(k,,) . (e^-ef +»/»?) . 

and q, v are the transferred momentum and frequency defined in (jdj). Following the principle to establish the contour 
time formalism in Section U and in the non-transient approximation^ (shown in lA II) . Li can be expressed in contour 
time formalism as 


Ti=- (T,5,VT(t')V2(t))o, 


(24) 


where is the time ordering operator defined in the Schwinger-Keldysh contour C and t G C+,t' G C_. The contour 
time evolution operator Sc is defined in contour C as 


Sc = Tce~^^-'^*^^'d)_ 

Trle-'^'^OAl 

In formula (1241) . the operators are defined in interaction representation by Hq and (A)o = — 

The resonant scattering probability r 2 is shown to follow 


(25) 




-pEi 


PiPf 


{-if 


2 ! 


E 


[^ 1 ^ 2 ] 




where 7r;(ti,t2) is given by 


T^i{h,t2)=Tt[Vn{h)Vi2{t2)]. (26) 

The /-dependent interactions Vn and Vi 2 (/ = 1,2) are defined by 

= Y e^(q,-)e—(27) 

^ ,/2eoVhujq^ 

Here qj,u!j,j = 1, 2 are all /-dependent defined in Table HI In the Schwinger-Keldysh contour time formalism, r 2 can 
be expressed in the non-transient approximation as 

^'^ = \Y f {TcSc'El{t'l,t2)'El{ti,t2))o, (28) 


where /i ,/2 G G C-. Following a similar procedure, the mixed scattering probability ri 2 can be expressed 

in the contour time formalism as 


ri 2 = Re 


*E / {TcScVl{t')Tri{ti,t2))o , 

, iJlt'u] 


(29) 


with /i,t 2 G C+,t' G C-. 
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TABLE I: Parameters q^j,ujj,j = 1,2 ofVij{t) in /-dependent. Also included are the /-dependent time factor Ti{t\,t 2 ). 


/ 

qi 

q2 

1^1 

6J2 


1 

-p/ 

Pi 

-UJf 



2 

Pi 

-P/ 

UJi 

-UJf 



B. Green’s functions in superconducting state 


In this article, we focus our study on the Raman responses of a superconductor in superconducting state. Thus we 
can simplify the electron Hamiltonian H = Ht + Hj in mean-field approximation as 

-^0 = 5I^k(^kr3 + AkTi)Tk, (30) 

k 

where d/k is the so-called Nambu spinor defined as iPk = (dkti band energy, and Ak is the paring 

potential of the Cooper pairs. Ti,i = 1,2,3 are the Pauli matrices. Introduce transformation iPk = /7(k)$k where 
d’k = (/kt) /-kj.) ) Hq can be diagonalized into the form as 

= $k, (31) 

k 


where ifk = vS^T+A^ is the diagonalized energy. Transformation matrix is defined as U(k) = Uk — ir^Vk, with the 
matrix elements Uk = 

Define the contour time Green’s function for the two-component operator $k as 

Ge(k;ti,/2) = -i(rc$k(H)$],(t2)). (32) 

For realistic calculation, we introduce the corresponding real-time Green’s function as 


(l + and Uk = sgn (Ak) (l- ^)- 


G(k;ti,/2) 


G^(k;G,/2) G<(k;G,/2)\ 

G>{k-,ti,t2) G^(k;/i,/2) j ’ 


(33) 


where an additional Schwinger-Keldysh index is introduced according to whether ti,t 2 belong to G+ or G_— . The 
four real-time Green’s functions in G(k;/i,t 2 ) are defined by 

G^(k;/i,t2) = ~*(d’k(G)‘hk(i2)), G^(k;ti,t2) = *(d’k(/2)‘hk(ii)), 

G'^(k;/i,t2) = 6>(ti - t2)G>(k;fi,t2) + ^*(^2 - H)G<(k;/i,t2), (34) 

G^(k;/i,t2) = 0(t2-ti)G>{'k;ti,t2)+0(ti-t2)G<{'k;ti,t2). 

Fourier transformed the real-time Green’s function by G(k;/) = ^fia;G(k, a;)e“*‘^‘ and G(k, w) = 

(i/G(k;/)e*“* with the real-time translational symmetry, it is readily shown that the real-time Green’s functions 
in mean-field superconducting state follow 


G<(k,..) 

G>(k,..) 

G^iKco) 

Gf(k,cc) 


ink^2TTS{uj — Ek) 0 \ 

0 m_k4.27r(5(w -|- Ek) j ’ 

—i (1 — rikf) 27r(5(a; — ifk) 0 

0 -i (1 - n-ki) 27r(5(a;-I-Flk) 

1—nkt J_ 7ikt 

uj — E\^-\-i6+ ' uj — Ei^ — i5+ 

0 

-(i-»kt) I -»kt 
bj — Ei^ — iS+ bj — Ek+iS+ 

0 


0 


1—n 


k4- 




a;+£^k+*<5"*‘ ~ io-\-Eiz—iS+ 

0 

-(l-n_kj.) I -rt-kj, 
0J+Ek—iS+ ' ui+Ek+iS+ 


( 35 ) 


















where rik-f = (/k-i./kt)o = = (/lki/-k 4 .)o = and d+ is a positive infinitesimal value. 

The reduced interactions Vij{t) and V 2 (t) defined in Section Til Al can now be re-expressed by the new two-component 
operator $k as 


= ^E‘^k+qWT^2(k,q)<l>kWe-*"‘, (36) 

k 


where q, u are the transferred momentum and frequency defined in Eq. © and qj,ujj{j = 1,2) are Z-dependent as 
given in Table HI Ti and T 2 are defined by 


Ti(k,q,) 


T 2 (k,q) 


1 

^2eoVhiOq- 




A(k, q) 


£k 

£Zk 


T 3 


Ak 

£Zk 


n 


(37) 


To obtain T 2 (k, q), we have considered the approximation q —>■ 0 in Raman scattering. Fig. |3] shows the schematic 
Feynman diagrams for the interaction vertices. 



FIG. 3: Feynman diagrams for the interaction vertices, 
factor T 2 (k, q). 



(a) for Vij with vertex factor Ti(k,qj) and (b) for V 2 with vertex 


C. Non-resonant scattering probability Fi 


Consider the non-resonant scattering probability Fi with the contour time formula (j24l) . In superconducting state, 
we consider a mean-field approximation where we neglect the roles of the collective low-energy excitations as we have 
discussed above. F 1 can be approximated in zero-th order mean-filed perturbation as 

r(°)=-/' {T,vUt')V2{t))o, (38) 

Jlt't] 

which is shown schematically by the Feynman diagram in Fig. 01 



FIG. 4: Non-resonant Fi in mean-field superconducting state. The solid line with arrow represents the Green’s functions (I35II 
and the vertex square represents the factor T 2 (k, q) in V 2 (t). 


Use Wick’s theorem to decompose the multi-particle correlation in F 


(0) 

1 


and then transform it into the real-time 
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formalism, the non-resonant scattering probability F® in mean-field approximation is given by 

q)G> (k + q; t' - t) 

k dti 

T2(k,q)G<(k;t-t')]e-*"(‘-‘'^ 

Ay: _ r 

= / ^^Tr[T^(k,q)G^(k-hq;a;-f i/)T2(k,q)G5^(ki;a;)]. (39) 

To obtain At in the last derivation, we have used the identity — v') = ), Substitute 

the Green’s functions Eq. (1551) into Eq. (15^ . F® is shown to be 



27rAt 

iV2 




c+-6{u + Ek - Ek+q) I 
+C— ^S{l' — Ek + Ek+q) \ 

+Cj^+5{v -\- Ek -|- Ek+q) 

-l-c— 5{v - Ek - Ek+q) ^ 


where c±± are defined by 


(40) 


E, 


c+- = 

'A 

Al 

C++ — (1 — n-k-q. 


2 (1 ~ ^k+qt) ^kti C —I- = (1 — U_k—qj.) n_k4., 

^1, 


K 


k 

A2 

If) nwt, c— = -^ (1 - nk+qt) n-ki- 


In formula (I40E c-|_and c_terms describe contribution from single-particle excitations and C++ and c_terms 

from Cooper pairs. At low temperature T Tc, since rik-f- = 0 and n_k 4 , = 1, c_|_= c_= C++ = 0 and c_=1. In 

this case, only Cooper pairs provide finite contribution to the non-resonant scattering, and thus 


■p(O) 
^ 1 


27rAi 

A2 





Ek — Ek+q). 


(41) 


It shows that there is a threshold frequency Vc beyond which the non-resonant scattering probability is finite. In a 
s-wave superconductor = 2A with A the superconducting gap. Our result Eq. (14111 is same to the previous one in 
the non-interacting limit and q —!> 0^^. 


D. Resonant scattering probability r 2 

In the mean-field superconducting state, F 2 can be approximated at zero-th order as 


do) 




Substitute Eq. (uni) of 7r/(fi,t2) into this formula, F^^^ is re-expressed as 


F^,°^ = 


E OM{k',q'})0;({k„q,})Ff\ 


47V2 


(42) 


(43) 


where qj,q), j = 1,2 are ^dependent given in TableUand Oi({kj,qj}) = Oi(ki, qi; k 2 , q 2 ) are defined by 


Oz({k„q,}) = Ti(ki,qi)Ti(k2,q2). 


In formula (H51) . F 2 * ^ is defined by 


■p(^^ ) _ 

i 2 — 


r;(t;,t')Tz(G,t2)(re$t (t')$k'+q'(4)‘i>l(i'i)$ki+q;(i'i) 

' + .1 ^ 


Crhf fd ' 


(44) 


‘^ki+qi (C)4’k^^_q^ (^2)4>k2 (^2))o, 


(45) 
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where ti,t 2 € (7+, € C- and the Z-dependent variables Ti{ti,t 2 ) = e *“^*1 are given in TablelH 

In the following, we will make all the Wick’s decompositions for the many-particle correlation in r 2 * \ They can 
be classified into three categories, the Rayleigh scattering, the fluorescence and the intrinsic energy-shift resonant 
Raman scattering. 


1. Rayleigh scattering 



FIG. 5: Feynman diagram for elastic Rayleigh scattering. 1,2 (T,2') on the vertices represent the time ti, t 2 (ti, t' 2 ) in the 
forward (backward) time contour C+ (C_). 


Rayleigh scattering is an elastic scattering with the incident and scattered photons having same frequency. The 
Feynman diagram for Rayleigh scattering is Fig. [5l which shows us that 


r 9 1 — :::: 


At 

27r 


Tr 


,T 


Gq (k'l, (k'l -I- -h 0j[) 


Tr [GJ(ki, vi)G ]^(ki + + wi)] 6r, 


(46) 


where dvidvi and 6r — <5k2-|-q2,ki<5q2,-qi^k(+q(,k(^q2,-q( • 

The Rayleigh scattering probability denoted by F^^J is shown to be 


rg = ^|/.r 

TT 


where the exchange symmetry between 1{V) = 1, 2 has been considered and Ii is defined by 

^ E / ^^OiTr [G^ (k, oj)G ^(k - p„ w - to,)] . 


(47) 


(48) 


Here Pi and ooi are the momentum and frequency of the incident photons and Oi = Oi(k, — pi;k— Pi,Pi). In the 
superconducting state, Ii follows 

(l-rikt)nk-pjt _ Wkt(l-"k-pit) 

uii+Ek-p,-Ek+iS+ uii+Ek-p,-Ek-iS+ 

I (1 kr )^—k+p^4- ri_k4,(l ^ —k+pj4-) 

— ^i — Ek—p, 

This formula shows clearly that only the single-particle excitations have dominant contribution to the elastic Rayleigh 
scattering in superconducting state. At low temperature T <C T^, the finite superconducting gap leads to /i = 0 and 
thus the Rayleigh scattering is strongly suppressed. 

It should be noted that since Ii = 0{N^) = 0(1), the Rayleigh scattering probability F^^j is also in order of 0(1). 

This is in contrast to the non-resonant scattering probability fJ°^ in Eq. (140]) or (|4T|) which is in order of 0(l/iV). If 
there is no special mechanism to suppress the Rayleigh scattering, it will be several orders of magnitude larger than 
the non-resonant scattering in contribution to the scattering cross section. 


k 


N 









11 



FIG. 6: Feynman diagram for fluorescence scattering. 


2. Fluorescence 


In a fluorescence process, there are sequent photon absorption and photon emission as schematically shown by the 
Feynman diagram in Fig. [5] Wick’s decomposition for the fluorescence scattering shows that 


_ 

^ 2.2 ~ 


2 X 


At 

27r 


f Tr [G^(ki+a;i)G^(ki,i/i)] 

J 1/1 I/O 


Tr [G^(k 2 +q 2 ,i ^2 +W2)G;5^(k2,i/2)] ^/, 


where <5/ = ^ki.k;<5k2.k^^qi,qi^q2.q2<^ti' s-nd are /-dependent as given in Table HI 

The fluorescence scattering probability denoted by F 2°2 is shown to be 

p(0) _ At (2) (2) 
i- 2 2 — - J1 1 J1 


2 ^ ldG.2, 

(2) 

where /-dependent integral t is defined by 

^1^3 I Ti (k, qj) I ^ Tr [Gs* (k -h q^, i/i -f ujj)G^ (k, vi)] . 

k 

( 2 ) 

In the mean-field superconducting state, !( ■ follows 


G = ^i:in(k,q. 


2 J nk'f'(l ^k+qj'f')^(^j T ./^k-t-qj li/k) 

1 Tn_ki,(l k—qji,)^(^_7 li/k+qj 


(49) 


(50) 


(51) 


(52) 


N ^ 1 -l-n_n fl — . I .-I-£/k) J 

It shows that only the single-particle excitations have contribution to the fluorescence scattering. At low temperature 
f2’) 

T <C Tc, II ■ = 0, thus the fluorescence response is largely suppressed in superconducting state. Moreover, the 
fluorescence scattering probability F^^^ is in order of 0(1) similar to the Rayleigh scattering 


3. Intrinsic resonant Raman scattering 

The intrinsic resonant Raman scattering processes are shown schematically in Fig. [7] and [5] They are all one-loop 
Feynman diagrams with four vertices of Vij and are classified into the two categories, ones shown in Fig. [7] where the 
times in contour branch G+ do not cross the times in G_ and the others shown in Fig. [5] with times cross. 

Denote the intrinsic resonant Raman scattering without time evolution cross by F^^^g. The four Wick’s decompo¬ 
sitions for r 2*^3 as shown schematically by the four Feynman diagrams in Fig. [7] have same contribution because the 
exchange of /(/') = 1, 2 is equivalent to the exchange of the time arguments. Therefore, we need only to consider one 
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FIG. 7: Feynman diagrams for intrinsic resonant Raman scattering without time cross. 


Feynman diagram such as Fig. [7] (a) with an additional factor 4. The Raman scattering probability F^^l 
be 


1 


p(o)__ 

2.3 4^2 




kill' 


where Of, = Of{ki + qi - qi,qi;ki - q 2 ,Pi - P/ - qi), Oi = Oz(ki,qi;ki - q 2 ,q 2 ), {q^, Wj, q', w'} , 
^-dependent given in Table ID F^^^ is given by 


ffj = -4 X ^ 

27r 


At 


Tr[G'^(ki + qi - 1 ^ 2 ) 0 ^(ki + qi,v'i) 


U-Uj 

'»< 


Go (ki,i/i)Go (ki - q2,J^2)](5ri/, 

where f , . = dv'iV 2 dvidv 2 and Sn, = 5 {v 2 — vi — loi + oj[)S{i''i — i^i — uji) 6 {iy 2 — + W 2 ). 

l/j Vj 

In the mean-field superconducting state, F^^^ follows 

^2.3^ = 87rAt[/34(5(w/ -UJ^+ ^^ki+qi - ^’ki-qs) 

T/3,2^(a^/ H^ki-l-qi 4" .^ki—q2)]; 

where are defined by 

1 - n-ki-rqi-q'jt , n-ki+qi-q'it 


-^3,1 = 


£'ki+qi + ^ki+qi—q'l + I?ki-|-qi + .E'ki+qi-q'^ 

1 — tlkit n-kit 


1^,2 — 


uji F^ki-t-qi 4“ I4kj^ iS~^ uj± F/ki-t-qi 4“ E]^-^ -t- id'^ 

(1 — riki-l-qit) n-ki-qatj 

1 ~ n._ki-qi+q'j^4, n_ki-qi+q'^4. 


-b 


4- Flki+qi ^'ki+qi—q'j 4- -f F/ki-l-qi .^'ki+qi-q'j 

1 n_ki4. 


4“ mki-t-qi F'k^ iS~^ uj\ “b F/ki-t-qi E'\^^ -b iS'^ 

(1 ~ n-ki-qij,) n_ki+q2^- 


is shown to 

(53) 

j = 1, 2 are 

(54) 

(55) 
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Similar to the Rayleigh and fluorescence scattering processes, only the single-particle excitations contribute to the 
non-time-cross intrinsic resonant Raman scattering. At low temperature T <C the finite superconducting gap 
strongly suppresses the scattering probability 




FIG. 8: Feynman diagrams for intrinsic resonant Raman scattering with time cross. 


From a similar derivation, the time-cross intrinsic resonant Raman scattering probability denoted by F^*^! is shown 
to follow 


■p(o) _ 

^ 2,4 — 


1 

IFp 


E OPO,Fg\ 

ki/r 


(56) 


where O*, = O*, (ki, q'^; ki — q 2 + Qi, Pi — P/ — Qi) and Oi(ki, qi; ki — q 2 -I- q'l, q 2 )- The two Wick’s decompositions 
shown by the two Feynman diagrams in Fig. [ 8 ] leads to 


-p(ii’) 
^ 2,4 


2 X 



Tr[Go>(ki-bqi,i/()Go<(ki,z.i) 


Go>(k2+q2,Z/')Go<(k2,i/2)]4„ki -q2-l-q'i^c5 


(57) 


where 5c = — vi — uji) 5 {iy 2 — vi — co[) 6 (iy 2 — I'l — uji + w^). 

In the mean-field superconducting state, F 2 4 ^ follows 


{ l4,l5{uJl -I- Gki — i?ki+qi)(5(w( -I- Gki — Aki+q'^) 

5{uJ2 - Wi - Aki -h Gki-q.+q;) ^ ^ 

+74,2<J(wi — Gki -|- Aki+qi)^(w( — Gki -I- Gki+q'^) 

5{uJ2 - wi -I- Aki - T^ki-qa-rq'i) 

where I 4 J are defined by 

-^4,1 = (1 ~ n.ki+qi'|-)nki^(l — nkj+q'j^-|-)nki-q2+q^t 

1^,2 = (1 ~ n._ki-qi^)n._ki4,(l — n_ki-q'j^4,)n_ki+q2-q'j4.- 

At low temperature T <C Tc, since nk-j- = 0 and n_k 4 , = 1, = 0. Therefore F® is strongly suppressed in 

superconducting state. Note that all the contributions from the intrinsic resonant Raman scattering to the cross 
section are in order of 0(1/Af), in contrast to that from the two-loop Rayleigh and fluorescence processes. 

In the general effective mass approximation for the electronic Raman scattering!^, the resonant and non-resonant 
responses are described by one uniform Raman charge density. In that approximation, the contributions from the 
resonant and non-resonant responses would be in proportion to each other. This is obviously in contrast to our results, 
where in superconducting state the non-resonant response has finite contribution from Cooper pairs while the resonant 
response has only contribution from single-particle excitations and is strongly suppressed by the superconducting gap. 


E. Mixed scattering probability ri 2 

The mixed scattering is a pure quantum effect, as it comes from the quantum interference of the resonant and 
non-resonant scattering processes. In the mean-field superconducting state, the mixed scattering probability denoted 
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FIG. 9: Feynman diagrams for resonant-non-resonant mixed scattering. 


by can be approximate by 


r® = Re 


. I J[t'U] 


(59) 


There are two Wick’s decompositions for as shown in Fig. [3 Because of the equivalence of the exchange of the 
time arguments and the exchange of I = 1,2, the two Feynman diagrams have same contribution to the scattering 
probability. Thus we only need to consider one Feynman diagram with an additional factor 2. F^'^^ is shown to be 


= Re 


iV2 


^ A* (k - q2, p* - Pf)Oi (k, qi; k - q2, q2 )r 


kZ 


(60) 


where (\j,j = 1,2 are ^-dependent given in Table U and is given by 


12 =2^^/ t^J^iTr[(^^r3 - ^ri ) G5’(k-f qi,i/i-hwi)Go (k,!/i) 


G^(k- q2,z/i - a;2)](5k',k-q2- 


(61) 


In the mean-field superconducting state, r [2 follows 


^(0 


where l 5 ,j,j = 1,2 are defined by 


= (dirAt) 


^k—c 


"[ ^k+qi ^k—q2 ) 


^k-q2 

“1“-^5,2-^k+qi ^k--q2 )] 5 


(62) 




h,l = (1 - «k+qit)»^k-q2t ( 
-^5,2 — (1 k—qij,)^—k+q2i 


~-E'k+qi+-E'k+^^'^ / 
_ l-n_k4. 

CJl+Sk+qj^ —Eiz—i 6 + 


-k4- 


\ t^i+-E’k-|-qi —Ek.-\-iS+ / 

It shows that only the single-particle excitations have contribution to the mixed scattering probability in supercon- 
ducting state. At low temperature T <C Tc, hj = 0, thus the mixed scattering probability r^ 2 ^ is strongly suppressed. 

It should be noted that r^ 2 ^ may be positive or negative in accord with the constructive or destructive interference. 
Moreover, it has magnitude in order of 0(1/N) similar to the resonant scattering with one-loop Feynman diagram. 


III. SUMMARY 

In the above sections, we have present a Schwinger-Keldysh perturbation formalism for the electronic Raman 
scattering. All the two-photon scattering processes can be well included within this contour time formalism and the 
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contributions from the resonant, non-resonant and mixed responses can be studied uniformly. As an example, we 
evaluate the Raman scattering cross section off an one-band superconductor. In the mean-field superconducting state, 
Cooper pairs contribute only to the non-resonant response. All the other responses from the Rayleigh scattering, the 
fluorescence, the intrinsic energy-shift resonant Raman scattering and the mixed scattering are dominated by the 
single-particle excitations and are strongly suppressed by the superconducting gap. The above formalism can be 
easily extended for the high-energy X-ray scattering when the inner core electrons are included to couple with the 
photon field^. A similar Schwinger-Keldysh formalism can then be established with a similar procedure. 


Appendix A: Review of Schwinger-Keldysh contour time formalism 

In this Appendix, we review the Schwinger-Keldysh contour time formalism which has been well established for 
non-equilibrium physics. This is a preliminary introduction for those who are not familiar with this formalism. More 
details can be found in Rammer’s textbook^. 


1. Contour time formalism 

Our task is to calculate the contour time correlation function defined by 

Oc = {Tc[AH{ti)BH{t2) ■ ■ ■ (Al) 

where Anit) = The Hamiltonian of the system is defined by 

H = Ho + Hi, (A2) 

where Hq is the quadratic part which can be treated exactly and Hj includes all the left such as the scattering 
potential and the inter-particle interaction, etc. 


(a) 

t 

; 

(b) tin 




(c) t in C 


t. 

I 


FIG. 10: Schematic illustration of the contour time representation of Anit). (a) for Anit) = S{ti,t)AHQ{t)S{t,ti), (b) and 
(c) for Auit) = rc[5cA_ffp (t)] with t € C+ and t £ C- respectively. 

Introduce the time evolution U matrix as and S matrix as S{t 2 ,ti) = 

then 

Anit) = U\t,ti)AHoit)Uit,ti) = S{ti,t)AHoit)S{t,ti), 
where Anoit) = In the contour time formalism, Anit) can be re-expressed as 


AH{t)=T4ScAHoit)], (A3) 

where Sc matrix is defined in the time contour C as 

Sc = Tce-i^-‘‘*^‘^*l (A4) 

Here Hi{t) = and the contour time ordering operator Tc is defined in contour C = C+ U C- 

with J^dt = dt. To obtain the contour representation (IA3I) . we have use the transitivity of the S matrix 
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Sit 3 ,ti) = S{t3,t2)S{t2,ti). For example, if t G C+, Anit) = S{ti,t)AHo{t)S{t,ti) = S{ti,tf)S{tf,t)AHoit)S{t,ti) = 
Tc[ScAHo{t)], and if f e C_, Anit) = S{ti,t)AHo{t)S{t,U) = S{ti,t)AHo{t)S{t,tf)S{tf,ti) = Tc[ScAHo{.t)]. The 
contour time representation (IA3I) is illustrated schematically in Fig. [TUI 

Following this principle, the contour time ordered correlation function Oc in (EU can be expressed as 

Oe = {T4ScA{h)B{t2) ■ ■ ■ Cit3)]), (A5) 

where the subscript Hq in the operator A, B and C has been ignored for clarity. 

c. 

c 

t,-ip 


FIG. 11: Three-branch contour (Kadanoff-Baym contour) Ci = C+ U C- U Ca, where Ca is the imaginary appendix contour 
defined as Ca ■ U ^ U — i/3- 


Now let us consider the thermal average. Since e = e with Ur = . Introduce an imaginary 

time axis, Ur can be expressed as 

Ur = Sc^ 

where the appendix contour Ca is defined as Ca ■ ti ^ ti — i/3 (shown in Fig. fTTll and the integral dt = dt. 

The correlation function Oc follows 


Tr 

e-^^°ScMScAiU)B{t2) ■ ■ ■ C(t3)] 

Tr 

e~d^oSc^ 



Tr 

e-P^°TcAScM)B{t2)---Cit3)] 

Tr 

e ^^°Sci 



(A6) 


where Tc^ is the time ordering operator in the contour Ci 


C+ U C- U Ca and matrix is defined by 


Sc. = Tc. 


j- dtH,(t) 


When we are not interested in the transient physics on the collision time scale, we can set ti —>■ —oo, then the 
contribution from the imaginary part of the contour Ci vanishes due to the thermal fluctuations^. The contour time 
correlation function can be approximated into 


Oc = {Tc[ScA{h)B{t2) ■ ■ ■ C(t3)])o. (A7) 

This is the non-transient approximation of the contour time correlation function. The time contour C in our study 
will be defined as the so-called Schwinger-Keldysh contour where ti —>■ —oo,tf —t +oo. (A)o is defined by (A)o = 
^Tr[e-P^oA]. 

When expand the Sc matrix in order of Hj and decompose the many-particle correlation function by Wick’s theorem, 
the perturbation corrections to Oc can be obtained order-by-order in principle. Wick’s theorem in the contour time 
formalism has been shown to have a same manner to the ground-state and the finite-temperature formalism. Define 
the single-electron Green’s function as 


Gc{l,2) = -i{Tcdidl), 


(A8) 
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where indices 1,2 involve the momentum, spin and temporal indices, etc. The zero-th order Green’s function is 
denoted by Go(l, 2). Wick’s theorem leads to all possible decompositions in an example as below: 

{Tcdid2d\d\)o 

= {Tcdid^^o{Tcd2d^^Q ± (T'c(iid3)o(T'cd2d4)o 

= -Go(l, 4)Go(2,3) T Go(l, 3)Go(2,4). (A9) 

where in last equation correspond to the bosonic and fermionic fields respectively. 


2. Real-time formalism 


The above formalism provides principle for the contour time perturbation theory. In realistic calculation, we will 
introduce the corresponding real-time formalism. In this formalism, the single-particle contour time ordered Green’s 
function Gc(l, 2) is transformed into a 2 x 2 matrix Green’s function G(I, 2), 


G(l,2) 


Gn(I,2) Gi2(1,2)\ 
G2i(1,2) G22(1,2) ) ' 


(AlO) 


where the subscribe indices n,m in G„rn(l,2) are the so-called Schwinger-Keldysh indices and are defined as n(m) = 
1, 2 according to ti{t 2 ) G G+ or G_. The real-time Green’s function has another familiar denotation 


G(l,2) 


G^(l,2) G<(I,2)\ 
G>(I,2) G^(l,2) ) ’ 


where the matrix element Green’s functions are defined by 

G>(1,2) = -i{did\), G<(1,2) = i(4di), 
G^(I,2) = -i{Ttdid\), G^(I,2) = -i{ftdid\). 


(All) 


(A12) 


It can be shown easily that 

G^(l, 2) = d(G - t2)G> (1,2) + e{t2 - ti)G<(l, 2), 
G^(l, 2) = e{ti - t2)G< (1,2) + 9{t2 - ti)G>(l, 2). 


The above formalism is defined for the fermionic filed. A similar formalism can be established for the bosonic field, 
where Bose-Einstein statistics should be introduced. Moreover the perturbation expansions in the real-time matrix 
formalism can be obtained one-to-one from the expansions in the contour time formalism. 
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